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FOREWORD 


This report comprises the analytical portion of a two-part 
vibration study of a pressurized torus shell. Part I of this study, 
which is contained in NASA CR-884, reports the experimental 
investigation and was prepared by Peter F. Jordan of the Martin 
Company, Baltimore, Maryland, under Contract No. NAS 1-6088. 
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VIBRATION STUDY OF A PRESSURIZED TORUS SHELL 


Part II — Development and Applications of Analysis 

by Atis A. Liepins 
Dynatech Corporation 


SUMMARY 


Previous work on the free vibrations of pressure prestressed toroidal 
shells is extended. Shells constructed of a number of toroidal segments, each 
with a constant thickness, and various boundary conditions are included. The 
analysis is then applied to simple models of a toroidal shell structure for which 
experimental results are published. The calculated fundamental frequency and 
mode shape of this shell agree well with the experimental results for three levels 
of internal pressure. 
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INTRODUCTION 


Several numerical analyses of the free vibrations of shells of revolution 
are reported in the literature. The most recent of these [ 1 ] includes the effects 
of prestress. Finite difference techniques are used in the numerical analysis. 

The formulation of the vibration problem in Reference [ 1 ] , however, is such 
that for thin toroidal shells the finite difference mesh converges slowly [ 2 ] , and 
for toroidal membranes it does not give meaningful results [ 3 ] . A different 
formulation of the problem for toroidal shells avoids these difficulties [2,3]. 

Previous analyses of the free vibrations of pressure prestressed toroidal 
shells [2,3] have dealt with unsupported shells with constant wall thickness. This 
report extends these analyses to toroidal shells with various support conditions 
and varying wall thickness. The supports are located at the inner and outer 
circumferences. The support conditions may be completely general. The shell 
may consist of a number of toroidal segments forming a complete torus. Each 
segment has a constant thickness. (Refer to Figure 1). All segments are made of 
the same material. 

The present analysis is applied to a toroidal shell structure for which 
experimental frequencies and mode shapes are published [ 4 ] . The structure 
consists of a toroidal shell stiffened at the outer circumference by a ring and at 
the inner circumference by a hub plate. Since the present analysis does not con- 
sider stiffening elements or additional masses, it is applicable to vibrations for 
which the stiffening elements are ineffective. Three simple models are investigated 
and the frequencies and mode shapes of the fundamental and first few overtones of 
each model are presented. 


2 



NOMENCLATURE 


\ 

h 

o 

P 


r 

u 

V 

w 


C, N, Q, S 


V T i> t 2* 


12 


D 


E 


E a’ E 0’ E q ! 0 


a 0 a6 

N , N , N . 
a 9 aO 

Q , Qn 

R 


S a’ s e 


a 


scale factor for mode shapes 

thickness of k-th toroidal segment 

reference thickness 

(1 - » 2 )pR/Eh k 

pressure 

( 1 - e cos a) / e 

meridional displacement (Figure 1) 

circumferential displacement (Figure 1) 

normal displacement (Figure 1) 

Eh /(I - v 2 ) 
o 

defined functions 

Eh 3 / 12 (1 - v 2 ) 

Young's modulus 

membrane strains 

difference spaces in k-th segment 

stress couples 

stress resultants 

transverse shear stress resultants 

radius of the generating circle of torus (Figure 1) 

membrane prestress forces 

meridional position angle (Figure 1) 

angle of k-th juncture 
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€ 



e 


K , K n , K . 

a 6 ctd 


X 


V 


p 

^6* ^ ctQ 
cv 



ratio of the two radii of the torus (Figure 1) 
h o /h k 

circumferential position angle (Figure 1) 

bending strains 
2 2 2 

(pR /Ee ) co, frequency parameter 
Poisson's ratio 
material density 

rotations of the normal to the shell 
circular frequency 
(I^/R) 2 / 12 

spacing between finite difference stations 


Matrices 


A, B, C, D, E, F, H 

P ’ X 1 ’ X 2 ’ Y 1’ Y 2 

I 

Z 

^ 1 *^ 2 * ^ 3 * ^4 


1x4 column matrices 
4x4 matrices 

Identity matrix 
1x4 column matrix 
4x4 matrices 


Indices 

i 

j 

k 

J 

K 

L 

n 


station in a segment 
station 
segment 
juncture station 
number of segments 
last station 

Fourier index of shell vibrations (corresponds to in 
Reference [ 4 ] ) 
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FUNDAMENTAL EQUATIONS 


The present analysis of the free vibrations of pressure prestressed 
toroidal shells is based upon a system of equations presented in Reference [ 5 ] . 
For the k-th toroidal segment these are: 

Equilibrium 


9 (rNJ + -§£?■ - N fl sin <* + rQ„ + 4" (1 + ~ L )Ss- <M„ 0 ) 


da ' a J 


a 2R 


r '80 'a0 J 


d E 


+ Ke> + (E a -E e )sina]+ rS^ ( -^ - 4> a ) 


+ pRr (f> a + ph^Rco ru = 0 


( 1 ) 


8N fl Q a 

- 7 T— + (rN ) + N sin o - Q cos a - [ ( 1 + — ) M . ] 

80 da ' O0 7 0 ^0 2R 8ck l ' r 7 o0 J 

+ r -^T- + 2E n sin a + d> cos a 1 + rS — - (E - + 0 .) 

0 L 80 o0 ^0 J a 8a ' a0 ^a0 / 


+ pRr <Pq + ph^R ct; rv = 


( 2 ) 


8 . SQ 0 . „ r dlp 6 


rN Q - N fl cos a - CrQ a ) - -gg- * S & \-gJ~ * <P a sin a - cos a] 

80 


+ rS ™ ( + El - pRr (E + E„) - ph. Rw 2 rw = 0 


(3) 


8 M 


aa < rM a> + ■ 

80 

8M » + 8 

(rM 

80 8 a 


Qf0' a 0 


(4) 


(5) 
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Strain- Di spl acement 


RE 


rRE 


9u 

a da 

~ dv 
e dd 


+ w 


+ u sin a - w cos a 


2rRE 


aO 
R K 


9v , 9u 

r *~r— + -z-z v sin a 


da d6 


d<p 
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a 


tRk 


da 

8 ( p 
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e 


de 


+ <P sin a 
a 


( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 


2rR/c 


aO 


r — + 


R <f> 


da 

9w 


a 


da 


dd " % sin « ’ ( r + cos a W ae ( n ) 

+ u (12) 


-rR 4> q = 


-2rR <p 


ad 


9 w 
de 

9u 

90 


+ v cos a 
9 


9 a 


(rv) 


(13) 

(14) 


Constitutive Relations 


Eh E 
k « 


Eh k E * = 


Eh, E . 
k a 0 


Eh k 
12 " a 

K 


K = 


12 0 

< K 

12 aO 


N - i>N fl 
a 0 


N, - ^ N 
0 a 


< 1+ N a* 


M - 
a 0 


M- - vM 
0 a 


(1 + v) M 


ad 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 
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The stress resultants S a and are known functions of pressure. 

They are determined from a separate analysis of the toroidal shell subjected 
to static internal pressure. An analysis based upon the linear membrane theory 
[ 6 ] gives : 


OL 


1 - 


€ COS OL 


~ pR 


1 - e cos a 


S 


Q 



( 21 ) 


REDUCTION TO SECOND ORDER DIFFERENTIAL EQUATIONS 


The solution of the fundamental equations (1 - 20) is started by separating 
the variables. Set 


~ - 

V *«• u 


r D i 

AN ,E , ~tP~ M , <P , 

o an an R an R an r an 

Ru ”1 
n 

_W M 9* V <*«• w _ 


D i D 

A N ,E q , ±-k q ,—g-Q , Rw 

L o On On R On R On R 2 H an nj 

N a0’ E a0’ ‘V ^ ae 



^o N «0n’ ^aOn’ ^On’ ^ aOn 


_ M aS ’ K a8’ %’ V _ 



IT^W T'W'P'Sn’ Rv n 



cos n 0 


sin n 9 
( 22 ) 


Next, define 


S = (sin a)/r 
C = (cos a)/r 
N = n/r 


(23) 


Then, use of equations (21), (22) and (23) in equations (1-20) yields (in nondimensional 
form) 

N: + S(N a - N > + NN„ a + 7,1 [Q„ + 4 ( 1 + C) M„J 


a 


ctQ 


k 'k L a 


, .. 2. 2^ 

+ (1 - v )e Xu = 0 


+ \ tl N < E a * - KJ + T S < E a ' E 0> + + i C >< E a- V + *«1 

(24) 

N a0 + 2SN a 9 - NN * ' r k”k f CQ « + T (1 + C)(M a* ‘ 

+ “k tl(-N E 9 + C* # )+ (l+i0)(E^+ *' a9 )+ <P e] 

+ (l-/)Av =0 (25) 
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N <* • CN 0 - Vk [(Qa + SQ a ) + NQ 9 ] + \ [t (N ^ + s<t> a ' C V 

+ (! + T C )( + E a ) - (E a + E fl )] - (1 - V 2 ) Aw = 0 (26) 

M i + S < M a " M 0> + NM afl ‘ Q a = 0 ( 27 ) 

M a« + 2 SM ae- NM e- % = 0 < 28 > 


E a 

— 

U* + W 

(29) 

E 0 

= 

Nv + Su - Cw 

(30) 

2E a 0 

= 

v T - Sv - Nu 

(31) 

K 

CL 

= 

4>' 

a 

(32) 

K e 

= 

S K + **9 

(33) 

2k 

a6 

= 

n - - N* a -(1 + C)* a9 

(34) 


= 

- w 1 + u 

(35) 

*9 

= 

Nw - Cv 

(36) 

2 <p n 
clB 

= 

Nu + v* + Sv 

(37) 

\ N a 

- 

E a + vE g 

(38) 

\ N e 

= 

E g + vE a 

(39) 

\ N ag 

= 

^-^ae 

(40) 

3„ 

V*a 

= 

K a + VK g 

(41) 

\ M 0 

s 

K e + VK a 

(42) 

3_ , 

Tj. M a 
'k aO 

= 

< 1 -") K a e 

(43) 
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In equations (24 - 43) and subsequent expressions, the subscript n has 
been dropped. Prime indicates differentiation with respect to o'. 

At this point the problem has been reduced to the simultaneous solution 
of equations (24 - 43). Next, we derive four simultaneous second order differ- 
ential equations for u, v, 0^, and Q^. The details of this derivation may be 
found in Reference [ 2 ] . The result is 


AZ" + BZ» + CZ = 0 


(44) 


where 


Z = 


a 


Q 


a 


and the elements of the A, B, and C matrices are given in the Appendix. 
Equations (44) are the governing equations for the k-th toroidal segment. The 
elements of the Z matrix are the basic variables of this formulation. 

The normal displacement is related to the basic variables by 
Qw = Tl u* + T 2 u + T 3 v- + T 4 v + T '« a + T 6 <P a - ^\<Q’ a + SQ a ) 

The Q and T’s are given in the Appendix. 


(45) 
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BOUNDARY AND JUNCTURE CONDITIONS 

An application of the principle of virtual work to equations (1-5) yields 
the boundary conditions. On an edge where a is constant prescribe: 


N + S E 
a a a 


or u — 0 


3 cos a 


N - (1 + ^ M + S (E +0 ) or v=0 

aO 2 R ' r ' aO or aO >0' 


9M . 

q + — _ s 4 > 

Rr d 0 a 


(46) 


or w ~ 0 


M 


a 


or 0 =0 

a 


Applying equations (21), (22), and (23) to the boundary conditions (46) yields: 


N « + ^r k k < 1 + i C)E 


a 


or u = 0 


N T. tk (1 + 3C)M n + — k (1 + ~ C)(E + 0 A 

aO 2 k k ' ' aS n, k ' 2 n aO ^a0 } 


Q + NM „ 
a ad 


kk 1 

v n + 


or v = 0 


or w = 0 


(47) 


r k^k 


M. 


a 


or 0 =0 

a 


where the subscript n has been dropped. 

The boundary conditions (47) may be written in a single matrix equation 


as: 


Vl + Y 2 y 2 = 0 


(48) 


where y^ and y^ are column matrices consisting of the nondimensional force and 
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displacement conditions respectively in equations (47); and are diagonal 
matrices that characterize the type of support condition. The elements of Y^ and 
Y^ matrices for two types of simple support, clamped support, and no support 
are given in the Appendix. 

It is of advantage to have the boundary equation (48) expressed entirely 
in terms of the basic variables Z. This is accomplished by direct substitution 
into (47). The result is 

y l = ^ 1 Z ' + ^2 Z 

(49) 

y 2 = V + *4 Z 


where the elements of the 0 matrices are given in the Appendix. 
With (49) the boundary equation (48) becomes 

t Y A + VsJ z ’ + t V 2 + VJ Z = 0 

At the juncture of two toroidal segments we require: 

— + 


(50) 


for equilibrium of forces 
for continuity of deformations 


y l = y l 
y 2 = y 2 + 


(51) 


where the plus and minus superscripts refer to the quantities just before and 
immediately after the juncture (see Reference [7]). 

With (49) the juncture equations (51) become 


<#>■ <Z')' + <*>~ Z' = </>* (Z + )' + <p + Z + 


(52) 


+ ’ 


<t> 3 ( Z ) + 0 4 Z = 0 3 (Z ) + 0 4 
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NUMERICAL ANALYSIS 


Since the geometry and prestress are symmetrical about a - 0 and 
a - 7r, we need to consider only one-half of the torus corresponding to the range 
0 ^ cl ^ 7T . Assume that one half of the torus is constructed of K segments, the 
k-th segment terminating at cl = /3 . Let each segment be subdivided by 1^ + 1 
equally spaced stations. Then the spacing between stations is 


A 1 = V 1 ! 

^k ' ^k-1 


*k 


*k 


k = 2, 3 . . .K 


and the position angle for the j-th station is 


(53) 


oi . - i A. 
J 1 


V p k-i + iA k 


k = 2, 3 . . . K 

i = 0, 1, 2 . . . ^ 


(54) 


The last station, a , is at a = 7 r, where 
L 


L = 


i + L i. 


k^l 


The derivatives of Z at all stations except those at junctures, Qfj , are 
approximated by the central difference formulas 


Z! 

J 





) 


zv 

] 



2Z. + 
J 


z 



(55) 
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With these formulas we obtain from equation (44) the set of difference equations 


D Z „ + E Z. + F. Z. = 0 (56) 

3 3+1 33 3 3-1 


where 


D. 

J 


E. 

J 


F. 

J 


A. + B. 

\ J J 





B. 

3 


(57) 


Equations (56) apply at all stations except those at junctures. 

In the juncture equations (52) the first derivative is approximated by a 
forward or backward difference formula 


< Z '>j = iT- (3Z j - 4Z J-1 + Z J-2> 
< Z+ >j = i^ (3Z J- 4Z J + l + Z J + 2> 


(58) 


Then for the J-th juncture, the juncture conditions in finite difference form are 


h 

2A 


£ 


+ „ + 


I < 3Z j - 4 Zj _ x + Zj_ 2 ) + f 2 Zj = - _ ( 3Zj - 4Z J+1 + Z J+2 ) + 0 2 Zj 


(59) 




2A 


L -< 3Z j' 4Z J-l + Z J-2> + ^4 Z J 


$ o 4- + + 

“+ (3Zj -4Z Ta1 + Z T+9 ) + 4> a Z 


J+l J+ 2' 4 J 

(60) 

At a = 0 the boundary equation in finite difference form is 


IS" *1 + Y 2 *2>< Z 1 - Z - 1> + (Y ? 0 3 + Z o " ° 


(61) 


and at a = tt 


+ Y 2 *2> < Z L + 1- Z L-1> + < Y >3 + Z L = 0 


( 62 ) 
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The superscripts o and ir denote the boundary matrices at or - 0 and 
a = 7r, respectively. 

The eigenvalues of the homogeneous equations (56) and (59) through (62) 
are found by trial and error using Potters method [ 8 ] , [ 2 ] . 


Z. = -P.Z.. , 
J J 3+ 1 


( 63 ) 


where 


-1 


P. = [E. - F. P. J D. 
3 L 3 3 3-l J 3 


(64) 


applies to all j except j = 0, j = J, and j = J+l. 


Write (56) for j = 0, and eliminate Z from (61). Then by comparison 


with (63) 


P o = [ 2A i( y!^ 3 + Y°0 4 ) + (Y°0 1 + ^0 2 )F' 1 E o ]- 1 (Y> 1 + Y° 2 4> 2 Hl * F o \) 


-1 


(65) 


Now, write equation (56) at j = J+l, equation (63) at j = J-2, and solve 
+ 

equation (60) for Z Then equation (59) becomes 
J 


B J Z J + 1 + Vj + F J z j-i = 0 


( 66 ) 


where 


D t = 


- 1 <** -X,X,^)(4I+ 


2A 


+ '^1 1 2^3 


= (+2 - X 1 X 2> + U- < *T - X lW 


F. = - 


<*1 -W3 )(4I+ P J-2> 


(67) 


and 


4> 

x i - - ^f + < 3i - D j+i F J+ i>]" 1 

+ ^ -i 

X 2 = *2 - ~ + < 3I - D J + lVl) 


The governing equation for the station after the juncture, where Zj is 
eliminated, is 
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0 


( 68 ) 


where 


+EZ +FZ 
J+l J+2 J J+l J J 


D J+ 1 D J+ l 


e j+i = e j+i ■ F j+i x i * 41 + D jii E J+ i> 


30 „ 4 > 

F t+1 = F X [<f> + ^ + -4- (41 +P T 9 )P T J 

J+l J+l 1 1 4 2A“ 2 A J-2' J-l j 


(69) 


Then Pj and Pj + ^ are found from equation (64) using (67) and (69). 

Equations (64) and (65) provide all the P's up to P 1 - Then write equation 

Xj - X 

(56) at j = L, equation (63) at j = L-l, and eliminate Z .. from (62). The 

Jj ' 1 

result is 


HZ t = 0 

Xj 


where 


(70) 


H - 2A K ( Y ^s + Y 2 0 4> ' < ^*1 + Y 2 *2> D L 1e 


- 1 . 


+ < Y 1*1 + V 2 >< 1 + d l p l )P l-i 


Since, in general, Z ^ 0, we must require that the determinant 

Xj 


V = H 


(71) 


vanish. Equation (71) is a frequency equation. 

When the vibration is axisymmetric (n = 0) and such that w is symmetric 
and u is antisymmetric about a = 7i> Z = 0. The frequency equation for this 

Xj 

case is obtained from equation (56) written at i = L-l with the help of equation (64). 
It is 


V = I E - F P 

1 L-l L-l L-2 


(72) 


For a natural frequency, is found from 


Z, = cH 
L, m m 


m = 1, 2, 3, 4 


(73) 
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where c is an arbitrary scale factor and the H are the cofactors of the first 

m 

row of H. The remaining Z’s are calculated in the reverse order with (64). The 
w component of the mode shape is calculated from equation (45) where the first 
derivative is approximated by the first of equations (55) at all stations except at 
the juncture and the station after the juncture. At these stations linear approxima- 
tions are used. 

The mode shapes are such that either v and w are symmetric, u is 
antisymmetric, or v and w are antisymmetric, and u is symmetric with 
respect to a = 0 and 7 r. These two groups of modes are denoted as symmetric 
and antisymmetric respectively. 

The foregoing computing procedure may be summarized as follows: 

1. Assume a value of X ; 

2. Calculate the elements of the A, B, C, D, E, F, and 
P matrices at all stations; 

3. Calculate the determinant; 

4. Repeat steps 1-3. From a plot of the determinant versus X 
determine a X (natural frequency) for which the determinant 

is zero; 

5. Calculate the mode shape for a natural frequency. 

The equations of this analysis are programmed in Fortran IV. The IBM 
7094 n produces a determinant in approximately 1. 5 seconds and a mode shape in 
3 seconds when seventy- five finite difference stations are used. Computation time 
increases linearly with the number of finite difference stations. 
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RESULTS 


The analysis developed in the preceding sections was applied to three 
simple models of a toroidal shell structure for which experimental frequencies 
and mode shapes are published [4j. The structure consists of a complete 
toroidal shell, made up of segments, each with a constant thickness; a plate, 
closing the hub; and a sheet metal ring at the outer circumference. A more 
detailed description of the model is given in Reference [ 4] . The average values 
of the structure parameters, taken from Reference [4] , are given in Table 1. 

Model I 

The weight of Model I equals the total weight of the toroidal structure 
including fittings. The weights of the hub plate with fittings and the outer ring are 
modeled as material added to the skin in the regions -0. 2^ a ^0. 2 and rr-0. 05^ a 
— ir + 0. 05, respectively. The width of these regions was chosen to keep the skin 
thicknesses reasonable. The wall thickness variation of Model I is shown in 
Figure 2. No attempt is made to model the elastic supports provided by the hub 
plate or the ring at the outer circumference. The model is supported in a free- 
free manner. 

The fundamental frequency and mode shape for 0, 6, 15, and 30 psi internal 

pressure was calculated using double precision arithmetic (16 digits). Convergence 

of the finite difference mesh was established for the frequency at 15 psi pressure 

with 181 difference spacings distributed as follows: 1^ = 11, l^= 18, I 3 = 28, 

I. = 119, I_ = 5. It was assumed that convergence for the frequency and mode 
4 0 

shape at other pressures was achieved with 181 spacings. Interestingly, the 
frequencies converged from below the final frequency. 

The overtones were calculated using single precision arithmetic (8 digits). 

A total of 75 difference spacings, 1^ = 1^ = 5, I^ = 10, I = 50, were used. 

The first overtone at 0 psi pressure was also calculated with 150 spacings in 
single precision and with 75 spacings in double precision. The frequency calculated 
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with 150 spacings was less than 1% lower than that calculated with 75 spacings in 
single precision. The frequencies calculated in single and double precision were 
identical to 3 significant digits. 

The calculated fundamental frequencies and mode shapes are compared 
with experimental results in Figures 3 and 4. At 15 psi pressure the calculated 
and experimental frequencies and mode shapes are in very good agreement. How- 
ever, the agreement deteriorates with decreasing pressure. The discrepancy 
between the calculated and measured results at 0 and 6 psi pressure may be due 
to an insufficient number of stations in the finite difference mesh. 

The calculated frequencies of the first five overtones of Model I are 
compared to the measured overtone frequencies in Table 2. There is wide 
disagreement between the calculated and measured frequencies. 

The mode shapes of the first five overtones of Model I for p = 0 are 
shown in Figure 5. The mode shapes for p = 6, 15, and 30 psi do not depart 
significantly from those for p = 0. 

Models II and HI 

Model II disregards the weight of the hub plate with fittings and the outer 
ring, as well as the elastic support provided by them. The wall thickness variation 
of this model is shown in Figure 6. Model II is supported in a free-free manner. 

Model in is an unsupported torus with a constant thickness of 0.0627 inches 
(h/R - 0. 00545). 

The fundamental frequencies and mode shapes of Model II were calculated 
with 150 finite difference spacings (1^ = = 1^ = 10, 1^ = 20, 1^ = 100). Again, 

the fundamental frequency converged from below the final frequency. All other 
frequencies and mode shapes of Models II and III were calculated with 75 spacings. 

The first five frequencies of Models n and HI for p = 0, 6, 15, and 30 psi 
and the measured results are compared in Table 3. The calculated mode shapes 
are displayed in Figure 7. Although the overtone frequencies of Models II and IH 
are closer to the experimental results than those of Model I, wide differences 
still exist. Furthermore, the large normal displacements at the hub of the com- 
puted symmetric modes are unrealistic because of the presence of the hub plate. 
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The effect of various boundary conditions on the first overtone of Model n 
was investigated. The following boundary conditions were considered: 


Hub (a = 0) 
simple support I 
simple support n 
clamped support 
elastic support 
free 
free 
free 
free 


Outer Circumference (a = ir ) 
free 
free 
free 
free 

simple support I 
simple support II 
clamped support 
elastic support 


In the elastic support the normal displacement and the transverse 
shear force were elastically related. The stiffness was taken to be that of the 
center plate under uniform in-plane load. The remaining three conditions in 
the elastic support were the appropriate ones of simple support I. 

The difference in frequency and mode shape for the elastic support 
and simple support I was found to be insignificant. This was expected as the 
in-plane stiffness of the center plate is high. 

The torus, in this mode, has identical frequencies and mode shapes for 
clamped support and simple support I. Computations verified this. The frequencies 
and mode shapes for the remaining four combinations of supports are shown in 
Figure 8. There is a small difference in frequency and mode shape when the torus 
has a simple support I, II and no support (see Table 3 and Figure 7 for frequencies 
and mode shapes of the free-free torus) at the outer circumference and the hub 
is not supported. However, there is a significant difference in frequency and 
mode shape when the hub is supported and the outer circumference is free. The 
importance of the conditions at the hub for this mode are further demonstrated by 
Figure 9. Here the normal displacements for Models I and II and the measured 
results are displayed. Model I and the measured results are in good agreement 
in the region outside the crown (90°^ a ^180° ). The two measured displacements 
inside the crown (O^o:— 90° ) fall between Model I and II results. 
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CONCLUSIONS 


Analysis and computer program are developed for free vibrations of 
complete toroidal shells with pressure. These shells may be constructed of 
several toroidal segments, each with a constant thickness. Supports may be 
specified at inner and outer circumferences . 

This analysis is applied to three simple models of a toroidal shell 
structure for which experimental results are published. Model I of this structure 
gives fundamental frequencies and mode shapes for three levels of internal pres- 
sure that are in good agreement with experimental results. Models more accurate 
than those analyzed here are necessary for the calculation of overtones. 

The effects of support conditions on the frequency and mode shape of the 
first overtone of Model II are briefly investigated. These calculations show that 
the support conditions at the outer circumference have a relatively small effect 
on the frequency and mode shape of this mode. However, careful attention should 
be given to the support conditions at the hub. 
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appendix 


The elements of the A, B, and C matrices are: 

A 11 = Q 2 [ 1 + k (1 + -tt C)] + QT 2 T ? 

A 12 = Q T 3 T 7 

A 13 = Q T 5 T 7 
A i4 = ' * S rQT 7 

A 21 " ^ T 1 T 9 
2 1 

A 22 ~ Q [ T (l-^) + k (l+Ic) + I(i- p ) r (i + 3C) 2 ] + qt t 

3 9 

A 23 = Q T 5 T 9 
A 24 = -^TQTg 

A 31 = 


A 32 

<© 

H 

co 

A - 

33 

QT, 

> 

> 

CO 

1) 

3 

-T? 

TQ 

A 41 = 

QTj 

T 

1 11 

> 

to 

1! 

QT 3 

T 

11 

> 

CO 

II 

Q 2 + 

« T 5 T 11 

A aa = 
44 

-. 3 

rQT u 
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B u = Q 2 S (1 + |k) + Q T ? (T^ + T 2 ) + T 1 T g 

B 12 = |Q 2 N [(1 + »)- j(l-») r (l +C)(1 +3C)] +QT 7 (T 3 ' + T 4 ) + T 3 T 8 
B 13 = QT 7 < T 5 + T 6> + T 5 T 8 

b 14 = -n 3 r< qst 7 + t 8 ) 

b 2 i = - |q 2 n[1 +v + |(1 - v) r (1 + 3 C) (1 - C) ] + Q Tg(T + Tg) + T 10 

B 22 = | Q 2 S [1-1- + k - T(1 +3C)(5 + 3C)]+ QT 9 (T 3 ' + T 4 )+T 3 T 10 

B 23 = TQ 2 N[^C+| (l-i'Ml + 3C)] + QT 9 (T g ' +Tg) +T 5 T 1q 

B 24 = -r, 3 r<QST 9+ T 10 ) 

B 31 = Q (T x ' + T 2 ) - Q' T x 

B 32 = « < T 3 ' + T 4> - <*' T 3 

B 33 = <* < T 5 + T 6> - T 5 

b 34 = -V 3 T(QS-Q') 

B 41 = Q T 11 < T 1 + T 2 > + T 1 T 12 

B 42 = - Q 2 N [ v C +|(l-v)(l+ 3C)] + QT 11 (T 3 ' + T 4 ) + T 3 T 12 

B 43 = « 2s + Q T 11< T 5 ,+ T 6> + T 5 T 12 

B 44 = -» 3 r<QST n+ T 12 ) 
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Q 2 [ (1 - v 2 ) eh + uC - S 2 - | (1 - i>)N 2 - |k (N 2 +S 2 ) - | (1 - (1 + C)2 N 2 ] 

+ Qt 7 t 2 ' + t 2 t 8 

Q 2 NS[-|(3-i’)-k+| (l-i/)T (1 + C)(1+ 3C)]+QT 7 T 4 ’ + T 4 T g 
Q 2 [kC + ~ (l-p)r(l + C)N 2 ]+ Q T 7 T 6 ' + T g Tg 
» 3 r[Q 2 -QT 7 (C-S 2 ) - S Tg] 

-Q 2 ns[|( 3-*) + k + ! (i-«)r(i+ O 2 ] + QT 9 T 2 ' + t 2 t 10 
Q 2 j(l - v 2 ) eh - | (1 -v) (C + S 2 ) - N 2 - | k (N 2 + S 2 + C 2 ) -kC - TN 2 C 2 

+ | (1- x)r(l + 3C)[S 2 (5+3C)-C(l+ 3Q]|+ Q T g T^' + T^g 
TQ 2 NS [C - | (1 - r)] + Q T g T g ' + T g T 10 

- r; 3 r [QT g (C-s 2 ) + st 10 ] 

- Q 2 + QT 2 ' - Q'T 2 

qt 4 '-q't 4 

Q 2 + QTg' - Q' Tg 

- V 3 r [Q(C - s 2 ) -Q'S] 

- | (1 - V) (1 +ON 2 Q 2 + QT n T 2 ' + T 2 T 12 

Q 2 NS [C +v (1 +C) + |<1 - v)(l + 30] + QT^' + T 4 
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to 1 1— 1 


where 


c 43 = Q 2 [vC-S 2 -|(X- ( .)N 2 ] + QT 11 T 6 - + T 6 T 12 

C 44 = +r [ QT X1 (C * s2) + ST 1 2 ]} 

T x = 1- pC + |kC-^(l-v)rN 2 (l-C)+|NT 0 [l + v+|(l-^)r(l+ 3C)(1 - C)] 

T 2 = S 0-C-k(l+ | C) + |(l-^)rN 2 ]+ NST 0 [|(3-|/) +k 

+ | (1- ^)T( 1 + C) 2 -|<l-^r (2 +30] 

T 3 = -|(1- i/)rNS-|-ST 0 [l-v +k + i(l-i/)r(l + 3C)(5 + 3C)] 

T 4 = -c - k(l+ C) + r [- N 2 c +| (l-i')S 2 -- (l-y)(l + 3C) c]j- 

- T 0 |(l - V 2 ) (1 - y)(C + S 2 )-N 2 -|k(N 2 + S 2 + C 2 ) -kC 

-TN 2 C 2 + |(1- v) r(l + 3C)[S 2 (5+ 3C) -C (1 + 3C)]j- 
T 5 = k (1 + ~ C) + TN 2 - TNT 0 [v C + | (1 - v) (1 + 3 C)] 

T 6 = S {| k + TN 2 + | rNT Q [3 (1- v) + (1 - 3v) C]| 

T ? = 1 - pC + k (1 + | C) + ~ (l-v) r<l +C) N 2 

T g = QS(1+C)[1+ “k - i (l - v) rN 2 ]-Q'T ? 

(1 - p) T NS (2 + 3 C) 

T 1q = QN [- v + C+ k(l + C) + TN 2 C + | (1 - i^) T (1 + 3C)(C-2S 2 ) 

- | (1 - p)T (1 - 3C)S 2 ] - Q 1 Tg 
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T u - | d + *)N 2 

T 12 = -(2N 2 SQ + Q'Tjj) 

T (1 - v) T (1 + 3 C) N 

X = 

0 | (1 - v) +k (1 + ±C) + | (1 - V) r (1 + 3 C) 2 

Q = (1 - v 2 ) e 2 \ - 1 - C 2 + 2vC - ~k (3 C + N 2 + C 2 ) - TN 2 [N 2 + (l-J^)C] 

+ T Q N j- v + C + k (1+C) +TN 2 C + | (1- y)r[(l + 3C)C-3<l+ C)S 2 ] | 

prime indicates differentiation with respect to a, and the subscript k has been dropped 
on r/, k, and T. The derivatives of N, S, and C are: 

N r = -NS 

S' = c - s 2 

C’ = - S (1 + C) 

Elements of ^matrix: 

T 

<t> u = + k (1 +|c) + [1- kC + k (1 + -|-C)] -J-} 

T 

* 12 = i [i- ,c + k< i + ic)]-^- 

T 

*13 = -jj-[ 1 -*C + k( 1+ §C)] 

0 14 = -r? 2 r[l-^C + k(l + |c)]^ 
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*21 

*22 

*23 

*24 

*31 

*32 

*33 

*34 

*41 

*42 

*43 

*44 


1_l; r(l + 3C)NS 1 


2 i ? - — Q 

[1+ ~r(l + 3C)(l+3C+4NS ~ (1 +|C) 

T, 


1_,; r(l+3C)NS 5 


2 7? 


Q 


-f (i-^)t? 2 r 2 (i + 30^ 


1 2 ^1 

-^3 d--)N^S-i 


4 T? 


ig (1 - v) N (1 + 3C + 4NS -^) 


~^S ^ 


(i - ») r 


n 2 s 

Q 


1 N 2 Tl 
— " N -Q 

n 


1 n 2 Ts 
"3 " N T 


V 


T 

2 ± 5 


X (1 + ‘'N -Q") 


- y 1ST T / Q 


Elements of matrix: 

T 

*11 = -jpf'S + [1-»C+ k(l+ I C)] — 
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T 

*12 = 4r N + [1-^C + k(X + 


*13 = ^[ 1 -f'C+ k(l + |o] 


Q 


014 = V r[l- *C + k <1 + i Q] -! 


Q 


0 


1 - v 


21 2 7] 


N [ 1 + J r (1 + 3C) (1 - C- 4S -|*] 


1 - V 


S [1+ i r ( 1 + 3 C) (1 + 3 C - 4 N )] 


22 2 1 ] 

23 = J YT r < 1+ 3qN(l + S^-) 


*24 =4 d-^Vd^oM 


0 


31 


1 (1 - v) N 2 (1 - C - 4S 


47? 


0 


32 


4 r\ 


V-( 1 - V ) NS (1 + 3C-4N -i-) 


0 


33 


- [(1 - v) N 2 (1 4 S-^) + ~r( 1 + |C)] 

71 ^ 


*34 = 1 + r d - 


. 1 2 T 2 
*41 = ^ -Q 
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(C - N 




-4 v ( S + N 2 
n 6 



) 



^ r 


n 2 s 

Q 


Nonzero Elements of tf> Matrix: 

O 

*31 = 

*32 - VQ 

*33 = V^ 

0 34 = - V 3 r/Q 

Nonzero Elements of <p 4 Matrix: 

*n = *22 = *43 = 1 

*31 = V Q 

*32 = T 4 /Q 
*33 = VQ 

0 34 =-*j 3 r s/q 

The subscript k has been dropped on n , k, and T. 
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Support condition matrices Y^ and Yg. Simple support I (v = 0): 


’ 0 

Y 2 " 

"l 

0 


1 

0 


1 

A. 




Simple support II (N^ = 0): 


Y 1 " 

r i 

0 

1 

Y 2 ’ 

"l 

0 


0 , 


1 


A 


_ |A-H_ 


Clamped support: 


No support: 


A 

Y = 

o 

”1a-i| 

1 A— 1| 

c 

A 

|A-l| 


A 

A 

i. J 


|A-1|. 


where A = 0 for symmetric modes 

A = 1 for antisymmetric modes 
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Example of elastic support: simple support I resting on an elastic foundation 
at the inner circumference 



k is spring constant 


The boundary matrices at a = 0 are: 


- 0 

Y 2 = 

‘ 1 

0 


1 

D 



- 0 3 |a-iJ 


1 

kR 


1 A— 1 1 , 


A 
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Table 1 


STRUCTURE PARAMETERS 


R 

€ 

E 

V 

PS 

Weight of outer ring 

Total weight of structure 
and fittings 


11. 5 in. 

0. 7326 

10. 7 x 10 6 psi 
0.32 

0.103 lb/in 3 
4.17 lbs 

76. 84 lbs 
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Table 2 


CALCULATED MODEL I AND MEASURED OVERTONE FREQUENCIES 


Calculated 

Mode 

Frequency, cps 

0 psi 
Pressure 

6 psi 

15 psi 

30 psi 

n = 1 
Symmetric 

471 

(178-190) 1 

476 

(212-224) 1 

482 

(185-201) 1 

490 

n = 2 

Antisymmetric 

550 

(252) 2S* 

577 

(301) 2S 

610 

(431) 2 

653 

n = 2 
Symmetric 

558 

(276) 2A* 

583 

(441) 2 or 3A 

615 

(444) 2 or 3 

653 

n = 2 
Symmetric 

562 

(391) 3S 

590 

(523) (3) 

625 
(521) 2 

671 

n = 2 

Anti s ymmetri c 

564 

(438) 2 

592 

(557) ? 

627 
(523) 3 

671 


The experimental frequencies (shown in parenthesis) are scaled from Figures 11a & 13 
of Reference [ 4 ] . 

* 2S - n = 2, Symmetric Mode ; 2A - n = 2, Antisymmetric Mode 
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Table 3 

CALCULATED MODEL H AND HI AND MEASURED FREQUENCIES 


Calculated 


n = 0 

Antisymmetric 
n = 1 


n = 2 
Symmetric 


n = 2 

Antisymmetric 


n = 3 
Symmetric 


n = 3 

Antisymmetric 


0 psi pressure 
II m Measured 


Frequency, cps 


6 psi pressure 15 psi pressure 

II IH Measured II HI Measured 



73.0 

48.8 

79.5 

93.1 

61. 2 

99.0 

115 

75.8 

123 

** 

178-190 ! 

** 

** 

212-224 

** 

** 

185-201 

** 


(!) j 



(1) 



(1) 


186 

252 

268 

219 

301 

296 

259 

431 

333 


(2S)* ** 



(2S) 



(2) 


199 

276 

274 

234 

441 

303 

276 

444 

343 


(2A)* 



(2 or 3A) 



(2 or 3) 


326 

391 

484 

372 

523 

526 

430 

521 

583 


(3S) 



(3) 



(2) 


325 

438 

490 

371 

557 

531 

429 

523 

587 


(2) 



(?) 



(3) 



* 2S - n = 2, Symmetric Mode; 2A - n = 2, Antisymmetric Mode 

** Calculated frequencies are higher than n = 3, Antisymmetric 










FIGURE I GEOMETRY AND NOTATION 
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AXIS OF 
ROTATION 



k 

JQ 

h 


k 

k 

1 

0.2 RAD 

1.800 

2 

30* 

0.160 

3 

58* 

0.100 

4 

17 - 0.05 RAD 

0.062 7 

5 

IT 

02687 


FIGURE 2 WALL THICKNESS VARIATION OF MODEL I 
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0 


4 5.4 


6 

15 

30 



48. 8 # 
61.2 ' 

75.8'* 


* AVERAGE OF TWO VALUES GIVEN IN FIGURE 9a OF 
REFERENCE [4] 

* * SEE DISCUSSION ON PAGE 10, REFERENCE [4] 

FIGURE 3 FUNDAMENTAL MODE OF MODEL I FOR n~0 
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VERTICAL DISPLACEMENT 






I ■ III! ■ IIB II II 


I 


l 


AXIS OF 
ROTATION 



n= I, SYMMETRIC 
471 CPS 




n * 2, ANTISYMMETRIC 
550 CPS 



n*2, SYMMETRIC 
562 CPS 


n*=2. SYMMETRIC 
558 CPS 



n» 2, ANTISYMMETRIC 
564 CPS 


FIGURE 5 OVERTONE MODE SHAPES OF MODEL I , p-0 
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XIS OF 
OTATION 



I 


AXIS OF 
ROTATION 


MODEL D 
MODEL IQ 



MODEL II FREOUENO 
GIVEN FIRST 


n = 0, ANTISYMMETRIC 
61.3, 73.0 CPS 



n * 2, SYMMETRIC 
242. 186 CPS 



n * 3. SYMMETRIC 
451, 326 CPS 



n = 2, ANTISYMMETRIC 
252, 198 CPS 



n «= 3. ANTISYMMETRIC 
458, 325 CPS 


FIGURE 7 MODE SHAPES OF MODELS II AND- HE. p-0 
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NASA- Langley, 1967 32 CR-885 


— MODEL I, 528 CPS 
MODELE, 246 CPS 



FIGURE 9: NORMAL DISPLACEMENT OF n= 2 SYMMETRIC MODE 
OF MODELS I SI, AND EXPERIMENTAL RESULTS 
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